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We report results on /b and semi-leptonic B decay form factors using NRQCD. We investigate 1/AI scaling 
behavior of decay amplitudes. For /b Effect of higher order relativistic correction terms are also studied. 



1. Introduction 

Weak matrix elements of B meson such as /b, 
Bb, and B ■n{p)li' form factors are impor- 
tant quantities for the determination of Cabbibo- 
Kobayashi-Maskawa matrix elements. However 
simulating the 6-quark with high precision is still 
a challenge in Lattice QCD, since the 6-quark 
mass in the lattice unit is large, amb 2-3, even 
in recent lattice calculations. One approach to 
deal with the heavy quark is to extrapolate the 
matrix elements for heavy-light meson around the 
charm quark mass region to the 6-quark mass as- 
suming 1/m scaling. It is, however, rather dif- 
ficult to control the systematic uncertainty in 
this approach, since the systematic error tends 
to become larger as increasing amg. An alter- 
native approach is to use an effective nonrela- 
tivistic action. In Table 0, we compare various 
features of NRQCDH, Fermilab[|, and ordinary 
Wilson/Clover actions as 6-quark action. The 
first two actions are the effective nonrelativistic 
actions. Their advantage is that 6-quark can be 
directly simulated. 

In this report, we present our study of the de- 
cay constant of B meson and B — > nW semi- 
leptonic decay form factors with NRQCD action. 
We study the mass dependence of these quan- 
tities by simulating heavy-light mesons over a 
wide range of the heavy quark mass. In section 
H, we study the 1/m dependence of the heavy- 
light decay constant using nonrelativistic action 
of ©(l/mg). The systematic errors due the trun- 
cation of higher order relativistic correction terms 
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are estimated. In section |^ we describe the first 
computation of the B ttW semi-leptonic decay 
form factors with NRQCD action of 0{l/mQ). 
Section ^ is devoted for discussion and future 
problems. 

2. B meson decay constant /b 

In the NRQCD approach it is very important 
to investigate the size of the systematic error aris- 
ing from the truncation of the action at a cer- 
tain order of l/mg. Earlier studies on fs by 
Davies et al.Q and HashimotoH] and the subse- 
quent work by NRQCD group g, where 6-quark 
is simulated with the NRQCD action including 
up to 0(1/toq) terms, showed that l/mp cor- 
rection of Jb from the static limit is significantly 
large. Thus the effect of the higher order cor- 
rection of 0(l/mQ) could be important. In this 
section, we compare the B meson decay constant 
obtained from the action including 0(l/m) terms 
only and that including 0(l/m^) terms entirely. 
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2.1. NRQCD action and field rotation 

We employ the following NRQCD action 



S = Q\t,x] 
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where n denotes the stabilization parameter. The 
coefhcients a are unity at tree level and should 
be determined by perturbatively matching the ac- 
tion to that in relativistic QCD in order to in- 
clude the 1-loop corrections. A and A'^^^ denote 
the symmetric lattice differentiation in spatial di- 
rections and Laplacian respectively and A*-"*) = 
^j(Ap'')2. B and E are generated from the stan- 
dard clover-leaf field strength. 

The original 4-component heavy quark spinor 
h is decomposed into two 2-component spinors 
Q and x after Foldy-Wouthuysen-Tani (FWT) 
transformation: 

Q{x) 



h{x) 
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where R is an inverse FWT transformation ma- 
trix which has 4x4 spin and 3x3 color indices. 



After discretization, at the tree level R is written 
as follows: 
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We apply the tadpole improvement to all link 
variables in the evolution equation and R by 
rescaling the link variables as ~^ U^/uq. 

We define two sets of action and current oper- 
ator {5H,R} as follows, 

setl = {5Hi,Ri} and scti = {SHn, i?//}(18) 
where 



5Hj = 5H^^^ and = 

1=1 

6 5 

SHii = ^^^'^ and Rii^Y. ^^'^ ■ 
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SHi and i?i include only 0{l/mQ) terms while 
SH2 and R2 keep entire ©(I/wq) terms and the 
leading relativistic correction to the dispersion re- 
lation, which is an ©(l/mg) term. The terms 
improving the discretization errors appearing in 
Ho and time evolution are also included. 

Using these two sets, we can realize two levels 
of accuracy of 0(1/toq) and 0{l/mQ). 

2.2. Simulation methods 

We have computed /s at /3 = 5.8 on 120 16"^ x 
32 lattices with periodic boundary condition in 
the spatial direction and Dirichlet boundary con- 
dition in the temporal direction. The inverse lat- 
tice spacing a^^ determined from nip is 1.714(63) 
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Figure 1. 1/Mp dependence of decay constant 
with the 0{l/mQ) action (open circles) and the 
©(l/mg) action (fiUed circles). 



GeV. For heavy quark, we use both mean-field im- 
proved 0(1/toq) and 0(l/m|) NRQCD action. 
We take six points for the heavy quark mass amg 
in a range 0.9-5.0 (1.5-8.5 GeV ). For light quark, 
we use Wilson action at k= 0.1600, 0.1585, and 
0.1570 (/ccrit — 0.16337) which correspond to mg- 
2ms. 

2.3. Results 

We show our results of fpy/mp in Figure ^ 
We find that the size of 0{\/m?p) correction is as 
small as about 3 % around the B meson region 
and about 15 % around the D meson region. 

In order to see how the correction terms in 
the rotation of operator (p^)-(|l6|) changes the re- 
sult, we show the contributions from each correc- 
tion term to fp\/Mp in Figure ||. We find that 
0{\/mQ) corrections are rather large. On the 
other hand, each of the three ©(l/mg) correc- 
tion is about 2 % around the B meson region. 
There is a cancellation among the three correc- 
tions, and the total eflFect is of 3%. Although the 
effect of the higher order corrections is naively 
expected to be very small, there is no guarantee 
whether this cancellation takes place at higher 
order. We, therefore, estimate an upper bound 
for the 0(l/mp) error to be of 6% at B meson 
region. 
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Figure 2. Corrections to the decay constant from 
each term of the current rotation. Circles rep- 



resent the correction of 1/toq term (13) with 
0{l/mQ) (open circles) and ©(l/mg)^ (filled cir- 
cles) actions. Squares, diamonds and triangles 
are corresponding to the corrections from 1/toq 
terms (|l^), (|l5|), and ([l^ ) respectively. 



We study l/mp dependence of fp^/mp by fit- 
ting the data with the following form 



fp^/rnp = [fp\/m^ 
for which we obtain 



rap nip 



p)°° = 0.308(20), 
ci = -0.87(11), 
C2 = 0.17(11). 

with the entire ©(l/mg) calculation. In physical 
units ci = -1.49(19)Gey,C2 = (0.71(23)GeV)2, 

To summarize, our analysis of Jb shows that 
©(l/mg) relativistic correction is of about 3 % 
and l/mp expansion from the static limit has 
a good behaviorl^. In order to obtain fp with 
higher precision, one has to control other system- 
atic errors such as perturbative and discretization 
errors. We have not included one- loop correc- 
tion for the renormalization constant, for which 
the calculation is underway. NRQCD group re- 
cently calculated the full one-loop renormaliza- 
tion factor for the decay constant with nonrela- 
tivistic heavy quark of ©(l/mg) and clover light 
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quarkQ. They find that the effect of operator 
dP^ which appear at one-loop level significantly 
reduces the decay constant. We are also planning 
to carry out the simulations at higher f3 values 
with 0(a)-improved Wilson light quark to remove 
0{a) error in near future. 

3. B meson semi-leptonic decay form fac- 
tors 

In this section we report our study of _B — s- TrlV 
semi-leptonic decay form factors. This is the 
first calculation with the NRQCD action. Ear- 
lier attempts to calculate the form factors were 
made by APE[|, UKQCDQ, WuppcrtalQ by 
extrapolating the results in the D meson mass 
region obtained with the clover action assuming 
heavy quark scaling law. Direct simulation would 
be certainly necessary as an alternative approach 
just as in the calculation of fs in order to in- 
vestigate the mass dependence and to obtain re- 
liable results. The status of form factors with 
Wilson/ Clover as well as Fermilab action is sum- 
marized in refs. |l2yi"^ ]. 

The semi-leptonic decay form factors and 
/° are defined as follows 



{7r{k)\V^\B{p)) 



p + k — q- 



'■B '"tt f0f^2\ 
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where = — fc^ and \B(p)) has a normaliza- 
tion 



{Bip)\B{p')) = i2nf2EBip)S^{p- p') 
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The virtual W boson mass takes a value in a re 
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gion < 0^ < where a'^ 
the rest frame of the initial B meson, pion is also 
almost at rest in the large q"^ region (g^ « 'Imax)^ 
where the W boson carries a large fraction of re- 
leased energy from the B meson. In the small q^ 
region (g^ Ri 0), on the other hand, the pion is 
strongly kicked and has a large spatial momen- 
tum (/Ctt ~ a few GeV/c). Lattice calculation 
is not reliably applicable for the small q^ region, 
since the discretization error of Olakj^) becomes 



unacceptably large. This leads to a fundamental 
restriction in the kinematical region where lattice 
calculation may offer a reliable result. 
The differential decay rate is given as 



dT{B° 



T+l- 



dq^ 
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where 

A('Z') = {ml +ml- qY - 4m|m2 . 
The decay rate vanishes at q^ — 



(24) 
because 

the phase space gets smaller. It is, thus, essen- 
tial to calculate f^{q^) in a q^ region where the 
experimental data will become available and the 
systematic error does not spoil the reliability, in 
order to determine \Vub\ model independently. 

Another important feature of the B meson 
semi-leptonic decay is the implication of the 
Heavy Quark Effective Theory (HQET)[|l|. In 
the heavy quark mass limit, it is more natural to 
normalize the heavy meson state as 

{M{p)\M{p')) = (2^)32 6'{p-p'){25) 
\ niM } 

instead of the covariant normalization ( p2[ ) . With 
this normalization, the large mass scale mu is 
removed from the theory and one can use the 
heavy quark expansion. The amplitude may be 
expanded as 

(^(fc)|t/jAf(p)) 

= (7r(fc)|F^|M(p))/V^ 

xfi + £l(!^ + £^ + ...) (26) 

where w is a velocity of the heavy meson and ci, 
C2, ... are functions oi v ■ k. It is worth to note 
that the heavy quark mass extrapolation and in- 
terpolation have to be done with v ■ k fixed. In 
the rest frame of the heavy meson, this condition 
implies fixed pion momentum, since v ■ k becomes 
\Jm'^ -I- p\. We propose to study the XjM de- 
pendence of the following quantities 

_ (7r(fc)|F4(0)|B(p)) 



V^{k,p) 



\J2Et^\J2Eb 



(27) 
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(28) 



which are natural generaHzation of fp^/ Mp for 
the the heavy-Ught decay constant. Indeed these 
quantities are almost raw number which one ob- 
tains in the lattice calculation as magnitudes 
of corresponding three-point functions, and then 
free from other ambiguities such as the choice of 
mass parameter of the heavy quark and the dis- 
cretization of the spatial momenta. 

3.1. NRQCD action and simulation pa- 
rameters 

The action we used for the semi-leptonic decay 
differs from that in the fs calculation. 
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Notation is the same as in section 2.1. The FWT 



transformation operator R is identical to that in 
eq. (|l^) Numerical simulation has been done on 
the same 120 gauge configurations as was used for 
the decay constant. For the light quarks we use 
the Wilson fermion at k = 0.1570. The follow- 
ing six sets of parameters for the heavy quark 
mass and the stabilization parameter; (mQ,n) 
= (5.0,1), (2.6,1), (2.1,1), (1.5,2), (1.2,2), and 
(0.9,2). TJiQ = 2.6 and 0.9 roughly correspond 
to b- and c-quark masses respectively. 

3.2. Extraction of three-point functions 

Matrix elements are extracted from the three- 
point correlation functions: 
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Figure 3. Effective mass plot for the pion with 
finite spatial momenta. 



Z^{k) Zb(p) 



2E4k) 2Eb{p) 

x(B(p)|T/;K(fc)) 

(for tB'>tv'> U) 



-tv) 



(32) 



\0b [tB , xbM (tv ,xv)Ol iU , 0) |0) 



where Ib , ty and indicate the location of the 
initial state B meson, heavy-light current and 
the final state pion respectively. We fix = 4 
and tv = 14, and tB is a variable. The ampli- 
tude {■K{k)\Vf,\B{p))/y/2E^y/2EB is obtained by 
dividing the above expression by the correspond- 
ing two-point functions. 

3.3. Effective masses 

In this section, we present our numerical re- 
sults. 

In order to see whether the contamination from 
the excited state is sufficiently small, we examine 
the effective mass plots of the two-point functions 
of pion and B meson, and three-point functions. 

Figure 3 shows the effective mass plot for the 
pion with finite spatial momenta up to \ak\ = 
2(27r/16). Correlation functions seem to reach 
to the ground state beyond t = 14, except for 
ak ~ (0, 0, 2) where statistical error becomes 
too large to extract reliable ground state en- 
ergy. We, therefore, use momentum values up to 
\ak\ = \/3(27r/16) in the following analysis. The 
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Figure 4. Dispersion relation for the tt and p 
mesons &t n — 0.1570. 
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Figure 5. The effective plots for the heavy-light 
two-point function at amq = 2.6 and k = 0.1570. 



maximum momentum value corresponds to ^ 1 
GeV/c in the physical unit. 

The momentum dependence of the pion energy 
should obey the dispersion relation 



(33) 



in the continuum limit. Then the deviation of the 
lattice dispersion relation from the continuum one 
is a good indicator of the discretization error. In 
figure ^ we plot the dispersion relation for tt and 
p mesons measured in our simulations. Solid lines 
represent the above expression (^3|), which does 
not fit the measured point, showing the effect of 
the 0(a) discretization error. 

A similar plot for the B meson is found in Fig- 
ure ^, where the maximum momentum is \ak\ — 
\/3(27r/16). The fitting interval is chosen to be 
16-24 where the contamination of excited state is 
negligible for each momentum values. 

Figures Hand show the effective mass plots for 
three point functions with temporal and spatial 
currents respectively. The location of the interpo- 
lating field of the B meson is varied to define the 
effective mass, then the energy should coincide 
with that obtained from the two-point function 
if the ground state is sufficiently isolated. We 



observe clear plateau in a wide range for t of 23- 
28 and the mass values consistent with that from 
the two-point functions, which are shown as solid 
lines in the figures. 
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Figure 6. Effective mass plot for the three point 
function with the fourth component of the vector 
current V4. 
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Figure 7. Effective mass plot for the three point 
function with the spatial component of the vector 
current Vj. 



Figure 8. I/'ti dependence of the matrix element 
V4 in eq.(|27|) at zero recoil. 



3.4. 1/m dependence 

We are now confident that the ground state is 
reliably extracted for both two-point and three- 
point functions, and present the results for the 
matrix elements. Since 1/m dependence of the 
matrix elements is the main issue in this study, 
we plot V4, and Vk , defined in and ( p8| ) respec- 
tively, which obeys the simple heavy mass scaling 
law ( 26 ) , in Figures || and||. We observe that 1/m 
dependence of matrix elements is rather small 
in contrast to the large 1/m correction for the 
heavy-light decay constant fp^mp. Although in- 
tuitive interpretation of this result is difficult, the 
smallness of the 1/m correction is a good news to 
obtain the form factor with high precision, be- 
cause the error in setting the fe-quark mass does 
not affect the prediction. And also this behavior 
is consistent with the previous works |p|, p^|JT^Jll[| , 
in which the heavy quark mass is much smaller 
than ours. 

The dependence of the form factors and 
/° is shown in Figure [l^. The heavy quark mass is 
roughly corresponding to the 5-quark. As we dis- 
cussed previously the accessible region is rather 
restricted. It is, however, interesting that already 
at this stage one is able to see the momentum de- 
pendence which could really be tested by looking 



at the momentum spectrum data in the future B 
factories. 
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Figure 9. 1/M dependence of the matrix element 
Vfc in eq. (Bq). Circles, diamonds and squares cor- 
respond to pion momenta k=(0, 0,1), (0,1,1) and 
(1,1,1) respectively. 
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is most reliable. 
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Figure 10. dependences of the form factors /° 



4. Discussion 

We have studied l/mq dependence of the 
heavy-light decay constant and semi-leptonic de- 
cay form factors with NRQCD action. We find 
that the error of truncating higher order relativis- 
tic correction term is as small as 6% for the decay 
constant. We also find that the semi-leptonic de- 
cay form factor for B nW has very small l/mg 
dependence, which is consistent with the previous 
results in the Wilson/Clover approach. 

To obtain the physical result for extracting Vut 
matrix elements, chiral limit for the light quark 
must be taken and calculation of the renormaliza- 
tion constant Zy at one-loop is required, which 
is now underway. 

One of the largest problem is that so far, due 
to low statistics and the discretization errors of 
0((afc)^), lattice calculation works only for rather 
small recoil region, where the statistics of the ex- 
perimental data is not high due to the phase space 
suppression. We are planning to carry out sim- 
ulations with much higher statistics, with larger 
(3 and with improved light quarks so that we can 
push up the accessible momentum region. It is 
also true that more data from CLEO, or future 
B factories is required for relatively small recoil 
region {k-^l ~ 1 GeV, where the lattice calculation 
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